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ABSTRACT 

We report the detection of GeV 7-ray emission from the molecular cloud complex that surrounds the su- 
pernova remnant (SNR) W44 using the Large Area Telescope (LAT) onboard Fermi. While the previously 
reported 7-ray emission from SNR W44 is likely to arise from the dense radio-emitting filaments within the 
remnant, the 7-ray emission that appears to come from the surrounding molecular cloud complex can be as- 
cribed to the cosmic rays (CRs) that have escaped from W44. The non-detection of synchrotron radio emission 
associated with the molecular cloud complex suggests the decay of tt^ mesons produced in hadronic collisions 
as the 7-ray emission mechanism. The total kinetic energy channeled into the escaping CRs is estimated to be 
Wesc ~ (0.3-3) X 10^^ erg, in broad agreement with the conjecture that SNRs are the main sources of Galactic 
CRs. 

Subject headings: acceleration of particles — cosmic rays — ISM: supernova remnants — radiation mecha- 
nisms: non-thermal 



1. INTRODUCTION 

Diffusive shock acceler ation (DSA) op erating at expand- 
ing shock waves (e.g., Mal kov & O'C Dru ry 2001) is widely 
accepted as the mechanism to convert the kinetic energy re- 
leased by supernova explosions into the energy of relativis- 
tic protons and nuclei (or cosmic rays) that obey a power- 
law type distribution. In DSA theory, cosmic rays (CRs) be- 
ing accelerated at shocks must be scattered by self-generated 
magnetic turbulence. Since the highest-energy CRs upstream 
in the shock precursor are prone to lack self-generated tur- 
bulence, they are expected to escape from the shock. The 
character of escaping CRs reflects complex interplay between 
CRs and magnetic turbulence, which is yet to be understood. 
DSA theory generally predicts that a substantial fraction of 
the shock energy is carried away by escaping CRs, and there- 
fore t hey should be treated as an integral part of the DSA pro- 
cess (lEUison & Bvkovll2011h . 

Depending on the amount of nuclear CRs released by a 
supernova remnant (SNR) and the diffusion coefficient in 
the interstellar medium, the molecular clouds illuminated by 
escaping CRs in the vicinity of SNRs are expected to be 
luminous in 7-ray s due to the enh anced tt^ -decay 7-ray s 
^haronian & Atovan 1996; Rodrig uez Marrero et al.l 120081 
lAharoniani2004i: iGabici et al.ii20Q9i). The TeV 7-ray emis sion 
seen in the vicinity of SNR W28 (lAharonian et al.ll2008bh can 
be regarded as a realization of this scenario. 
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Recent detections of intense GeV 7-ray em ission from 
SNRs interacting with molecular clouds (e.g.. lAbdo et al.l 
IMM IMoa bllcl : ICastro & SlandlMol: iGiuhani et al.ll2()TTh 
provide an interesting opportunity to study how CRs escape 
from SNRs and propagate in the interstellar medium. How- 
ever, it has not been clear whether the 7-ra y emi ssion is pro- 
duced by escaping CRs or by trapped CRs dUchiyama 20TTh . 
since shocked molecular clouds inside a SNR also could be 
the sites of efficient 7-ray production (iBykov et al.l 120001: 
lUchivama et"an[20TQh . 

In this Letter, we report the Fermi Large Area Telescope 
(LAT) observations of the region around SNR W44. We 
demonstrate that the 7-ray emission distributed on a spatial 
scale much larger than the remnant is attributable to a CR halo 
around SNR W44. 

2. OBSERVATION AND ANALYSIS 

The LAT onboard Fermi is a pair-conversion 7-ray detec- 
tor that cov ers a very wide ran^ e of energy from 20 MeV to 
>300 GeV (lAtwood et al.ll2009l) . The LAT tracks the electron 
and positron resulting from pair conversion of an incident 7- 
ray in thin high-Z foils, and measures the energy deposition 
due to the subsequent electromagnetic shower that develops in 
the calorimeter. The point-spread function (PSF) varies with 
photon energy and improves at higher energies. The 68% con- 
tainment radius is smaller than 0.°5 above ~ 2 GeV. 

We use the 7-ray data acquired from 2008 August 4 to 201 1 
September 6. The Source class events are analyzed using the 
instrument response functions (IRFs) P7SOURGE_V6. A 
cut on earth zenith angles greater than 100° is applied to re- 
duce the residual 7-rays from cosmic-ray interactions in the 
upper atmosphere. We analyze photons with energy > 2 GeV, 
where the PSF is sharp enough to disentangle multiple spatial 
components. In-depth analysis below 2 GeV will be published 
elsewhere. 

Spectral and spatial parameter estimation is done by maxi- 
mizin g the likelihood of the source model (e.g., M attox et aH 
119961) using "binned" gtlike of the Fermi Science Tools. The 
region used for the likelihood analysis is 10° x 10°, centered 
on W44. The 7-ray source model includes point sources listed 
in the second Fermi LAT catalog (2FGL sources: .Nolan et al.l 
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Fig. 1 . — (a) Fermi LAT 7-ray count map for 2-100 GeV around SNR W44 in units of counts per pixel (0° 1 x 0° 1) in celestial coordinates 
(J2000). Gaussian smoothing with a kernel a = 0°3 is applied to the count maps. Green contours represent a 10 GHz radio map of SNR W44 
JHanda et al.ll 19871) . 2FGL sources included in the maximum likelihood model are shown as crosses, while those removed from the model are 
indicated by diamonds, (b) The difference between the count map in (a) and the best-fit (maximum likelihood) model consisting of the Galactic 
diffuse emission, the isotropic model, 2FGL sources (crosses), and SNR W44 represented by the radio map. Excess 7-rays in the vicinity of 
W44 are referred to as SRC-1 and SRC-2. 

sky model. Significant excess 7-rays are seen in the vicinity 
of W44; the features are referred to here as SRC-1 and SRC- 
2. The statistical significance is found to be ~ 9a for SRC-1 
and - lOcr for SRC-2. 

The residual count map depends weakly on the choice of the 
spatial template that describes 7-rays from W44. Our simula- 
tions using gtobssim verified that SRC-1 and SRC-2 are not 
caused by photons leaking from W44 due to the PSF of the 
LAT. Also we checked the robustness of the results by select- 
ing only the front-converted events. 



l2Q12h , Galactic interstellar diffuse emission, and an isotropic 
component (extragalactic and residual particle background). 
The Galactic diffuse emission is modeled using the standard 
ring-hybrid model, gal_2yearp7v6_v0.fits, with its normal- 
ization being left free. We use a tabulated spectrum written in 
iso_p7v6source.txt as the isotropic diffuse emission. The 
LAT data, analysis software, and diffuse models are made 
publicly available through the Fermi Science Support Centei0. 

Figure [TJa) shows a 2-100 GeV count map in the vicin- 
ity of SNR W44, where crosses and diamonds indicate the 
positions of 2FGL sources. In addition to W44, five 2FGL 
sources are distributed within l.°5 from W44. One of them, 
2FGL J1857.6+0211, coincides with PSR B1855+02 and also 
with S NR G35.6-0.4 B that has recently been re-identified as 
a SNR (lGreenll2009h . The other nearby 2FGL sources (dia- 
monds) do not have clear counterparts in other wavelengths, 
and they are excluded from the source model to investigate 
the surroundings of W44. 

W e employ a synchro tron radio map of SNR W44 taken 
from iHandaet aH (Il987h to model the spatial distribution of 
the 7-ray emission from W44, given that the synchrotron and 
7-ray emission from W44 are expected to be co-spatial (see 
^4.1t . The 7-ray spectrum is assumed to obey a power law. 

3. RESULTS 

The likelihood analysis is performed using the source 
model described above. For point sources, we use the spectral 
models adopted for each source in the 2FGL catalog anal- 
ysis. Spectral normalizations of point sources located < 3° 
from W44 are allowed to vary in the likelihood fit, while the 
spectral parameters of the other field sources are fixed using 
the 2FGL catalog. The normalization and photon index of 
W44 are left free; a photon index of rw44 = 2. 94 ± 0.07 is 
obtain ed in agreement with our previous work (lAbdo et al.l 
l2010ab . Figure [lib) shows a residual count map, where the 
observed count map in 2-100 GeV is subtracted by the best-fit 

^ http://fermi.gsfc.nasa.gov/ssc/ 

10 PSR B1855+02 is located near the center of G35.6-0.4. At the south- 
ern border of G35.6-0.4, there is a TeV 7-ray source HESS J1858-h020 
(Aharonian et al. 2008a), toward which one or more molecular clouds have 
been found (Paron & Giacani 2010). Discussion of HESS J1858-h020 can be 
found in Torres et al. (2011). 
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Fig. 2. — (a) Fermi LAT spectrum of SRC-1 (red points) along 
with the LAT spectra of SNR W4 4 from this work (black points) and 
previous one dAbdo et al.ll2010al , gray points). The 7 -ray spectrum 
reported by AGILE is also shown (iGiuliani et al.l201 iT open circles). 
Systematic errors are added in quadrature to the errors of the SRC-1 
spectrum. The model curves describe the emission from SNR W44 
(see ^4.1K consisting of tt^ -decay 7-rays (dashed curves) and rela- 
tivistic bremsstrahlung (dotted curves), (b) Same as (a) but the LAT 
spectrum of SRC-2 is shown instead of SRC-1. 

We perform spectral analysis of SRC- 1/2 by modeling each 
source as a disk with a 0.°4 radius (see Fig. [TJ. The re- 
sulting 7-ray spectra are plotted in Figure |2] along with the 
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spectrum obtained for SNR W44. Adding the SRC- 1/2 disks 
to the source model does not significantly affect the W44 
spectrum. The power-law photon index is found to be F = 
2.56 ± 0.23sta ± 0.2sys and T = 2. 85 ± 0.23sta ± 0.2sys for SRC- 
1 and SRC-2, respectively. The systematic errors are evalu- 
ated from different choices of the sky models describing W44 
and SRC- 1/2 and from the uncertainties of the effective area. 
The imperfection of the diffuse emission model and its pos- 
sible impact on the results are discussed below. We tested a 
smoothly broken power-law and exponentially cutoff power- 
law for SRC- 1/2 but found that the spectral fits do not signifi- 
cantly improve. 




Fig. 3. — Fermi LAT residual count map highlighting the 7-ray 
emission from the surroundings of SNR W44. Magenta contours 
present the synchrotron radio map of SNR W44. Green contours 
show CO (7=1^0) emission integrated over velocity from 30 to 65 
km s~^ with respect to the Local Standard of Rest (iDame et al .12001) , 
tracing the molecular cloud complex that surrounds SNR W44. The 
contours start from 20 K km s~^ with an interval of 10 K km s~^ 
Some molecular clouds not associated with W44 are also seen along 
the Galactic plane (dashed line) in the CO map. 



4. DISCUSSION 

We have discovered GeV 7-ray sources on the periphery of 
SNR W44. It has long been known that a complex of giant 
molecular clouds (GMCs) surrounds SNR W44; the spatial 
extent is as large as 100 p c and the to tal mass of the co mplex 
amounts to -- 1 x lO^M© (iDame et al.n 1986 : Seta et al.|[T998) . 
In Figure[3j the 7-ray emission from the su rroundings of W44 
is compared with a CO (/ = 1 ^ 0) map (iDame et al.ll2001h 
integrated over a velocity range of 30 to 65 km s~^ appropriate 
for the GMC complex. The regions of excess 7-rays overlap 
with the surrounding GMC complex. 

The 7-ray emission in the vicinity of W44 can be ascribed 
formally to possible imperfection of the maps of gas column 
densities used in the model of the Galactic interstellar diffuse 
emission, or to a local enhancement of CR density. The for- 
mer implies that the mass in the 7-ray-emitting region around 
W44 is underestimated by a large factor (^5), or it requires 
the presence of unknown background clouds with a huge mass 
of > 1O^M0 . Therefore an overabundance of CRs in the vicin- 
ity of W44 offers a more sensible explanation; we present a 
model in which the GMC complex is illuminated by CRs that 



were produced in SNR W44 and escaped from it. We first 
model the 7-ray emission from W44 itself, and then proceed 
with modeling the 7-ray emission from the surroundings. 

We adopt the follow ing parameters to describe SNR W44 
(lUchiyama et al.ll2QlQl and references therein): (i) J = 2.9 kpc 
as the distance to W44 b ased on the firm association with the 
surrounding GMCs (e.g. JSeta et anil998h : (ii) R = 12.5 pc as 
the radius (corresponding to the angular radius of ^ = 14.8'); 
(iii) the kinetic energy released by the supernova Esn = 2x 
^q51 gj,g. ^j^^ ^j^g ejecta mass Mqj = 2Mq\ and (v) the remnant 
ageQ ^age = 10000 yr. Assuming evolution in the uniform in- 
tercloud medium, these parameters imply that the intercloud 
medium has hydrogen density of ^ ::^ 2 cm~^ and the Sedov- 
Taylor phase started around ^ = ^st — 129 yr when the radius 
was rsT — 1.9 pc. 

4.1. Gamma-ray Production Inside SNR W44 

A high-resolution radio continuum map of SNR W44 is 
dominated by filamenta ry structures of synchrotron radiation 
(ICastelletti et al.l l2Q07h . The radio emission is thought to 
arise from radiatively-compressed gas behind fast dissocia- 
tive shocks dr iven into molecula r clouds that are engulfed by 
the blastwave (iReach et al .120051) . Assuming typical magnetic 
fields of molecular clouds, the GeV 7-ray flux relative to the 
radio flux is expected to be high enough to account for the 7- 
ray emission from SN R W44, irrespective of the origin of the 
high-energy particles (lUchiyama et al.ll2QTQl) . 




Fig. 4. — Modeling of the 7-ray emission from the molecular cloud 
complex that surrounds W44. Data points are from Fig.|2] but SRC- 
1 and SRC-2 are co-added. The SRC- 1-1-2 spectrum (red points with 
statistical errors) is attributable to the vr^-decay 7-rays from the cloud 
complex illuminated by the CRs that have escaped from W44 (blue 
curves). Three cases of diffuse coefficient, D2S = 0.1, 1, 3, are shown. 
A gray curve indicates the 7-ray spectrum produced by the sea of 
GCRs in the same CR-illuminated clouds. 

The dense radio-emitting filaments are indeed the most 
probable sites of the dominant 7-ray production i n SNR W44, 
Riven the estimated mass of M^h = 5 x lO^M© (iReach et al.l 
120051) which is ~ 9 times larger than the swept-up intercloud 
mass. Assuming a pre-shock cloud density of no = 200 cm~^ 

^^ The SNR age is comparable to the spin-down a ge of the pulsa r 
B 1853+01 associated with SNR W44, r^d = 20000 yr (Wolszczan et al.ll991l) . 
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(iReach et alJ 120051) and a pre- shock magnetic field of ^o = 
30 fiG, we can estimate the compressed gas density and mag- 
netic field in the filaments as rim ^7 x 10^ cm~^ and Bm — 0.^ 
mG following the prescription in lUchiyamaet al. (2010). 

For simplicity we describe the energy distributions of CR 
electrons and protons in the filaments as a cut-off power law 
in momentum: rie^pip) = K.pP'^'^^ ^W(~P / Pc), where the in- 
dex is ch osen to match the radio spectral index of a ::^ 0.37 
(fCastellet ti et al.ll2007h . The ratio of radio and 7-ray fluxes 
yields ke/kp = 0.05. The spectral break in the LAT spec- 
trum is reproduced by pc = 10 GeVc"^ As shown in Fig. [21 
the spectrum below a few GeV is dominated by the decays 
of TT^ -mesons produced in the dense filament^ while the 
falling part of the LAT spectrum is contributed largely by the 
electron bremsstrahlung. The energy density of the CR pro- 
tons amounts to Up c^ 4 x 10^ (Msh/5 x 10^ M©)"^ eV cm"^, 
which could be explained by reacceleration of Galactic CRs 
(GCRs) pre-existing in a molecular cloud (lUchiyama et al.l 
I2OIO). However, we do not specify the dominant source of 
the 7-ray-emitting particles, since freshly accelerated CRs 
could also enter the filaments diffusively from the intercloud 
medium. Provided that Up corresponds to the mean CR den- 
sity in the shell of the remnant, the total kinetic energy of CRs 
is Wp - 0.4 X 10^^ (Msh/5 X 10^ Mor^ erg. 

4.2. Gamma- rays from the Surroundings ofW44 

The CR distributions in the surroundings of SNR W44 
should be determined by (1) how CRs are released into the 
ambient medium, and (2) the di ffusion coeffi cient of the in- 
terstellar medium, Dism(p) (e.g. JGabici et al.l l2Q09). We con- 
sider only CR protons since leptonic emissions are unimpor- 
tant for an e-p ratio of ~ 0.01. Also, the non-detection of 
synchrotron radio emission from the GMC complex indicates 
that electron bremsstrahlung is not the main 7-ray production 
mechanism. 

Let us assume that CRs with a momentum p can escape 
from the surface of a SNR at time t = tQ^dp) when the SNR 
radius becomes Rqsc(p), and that tQ^dp) has a power-law form: 



lowing solution of the diffusion equation (lOhiraet al.ll2011h 



tcscip) = tsT 



-1/x 



(1) 



where we a d opt p ^ax = 10^^ eVc | and y = 3, following 
iGabici et al.l (l2009h and lOhira et al.l (1201 Ih . This impHes 
that SNR W44 is currently releasing CRs with /? = /^o — 
2 GeVc"^ A simple Sedov-Taylor evolution gives Rqsc(p) = 

rsjip/Pma.r^^'''- 

The runaway CR spectrum integrated over SNR ex- 
pansion is expected to be of the form Nq^c(p) ^ P~^ 
(iPtuskin & Zirakashvilill2005|) . However, depending on the 
time history of acceleration e fficiency and maximum energy, 
it cou ld be different from p~^ (lOhira et al.l2010l:ICaprioli et~an 
l2010b . We parameterize the total spectrum of CRs injected 
into the interstellar space as Nqsc(p) = hscP'^ ^^p(-p / Pmax)- 

The distribution function of the runaway CRs at time t at 
a distance r from the SNR center, n(p^r^t), is described by 
a well-known diffusion equation, which can be solved using 
the method developed by lAtoyan et al.l (Il995|) . We use the fol- 



^^ The AGILE spectral data dGiuliani et al.ll20lll) are not taken into ac- 
count in the model. They will be discussed elsewhere in light of a low-energy 
spectrum measured with Fermi. 



n(p,r,t) = 



Ncsc(p) 

47ry^RdR,scr L' 



-(r-/?esc) /Rj 



-(r+Rescf/R] 



where 



Rd(pj) = 2^/DlSM(p)[t-t^sc(p)l 



(2) 



(3) 



The diffusion coefficient of the interstellar medium is often 
parameterized as 



Dism(p) = 10 



28 



/)28 ( 



10 GeV c 



cm^ s 



(4) 



with constants of Z) 28 ^ 1 and S 



06 based on the GCR 
The diffusion co- 



propagation model (iPtuskin et al.l l200i 
efficient in the close vicinity of a SNR may be different from 
the Galactic average, for exam ple, because of A lfven waves 
generated by CRs themselves (IFujita et aH 120101) . We allow 
D2S to vary, but fix the index to be ^ = 0.6 for simplicity. 
Given that the 7-ray emission by the escaping CRs is visi- 
ble against SNR W44 in the LAT image, we expect the size 
of the CR halo, Rcr(p) = ^j(/^,^age) + ^esc(/^), is much larger 
than the size of the remnant R = 12.5 pc. Since we find 
7?cr(100 GeVc"^) ^ 2R for Z)28 = 0.1, we set Z)28 > 0.1. Note 
however that the constraint depends on the exact choice of 
^age; ^28 > 0.05 is obtained for ^age = 20, 000 yr. 

To calculate the 7-ray emission from the surroundings of 
W44, we also need to specify the mass distribution. The total 
mass of the cloud complex illuminated by the runaway CRs 
is estimated as Muc ~ 5 x l O^M^;^ by adding up six molec- 
ular clouds (Clouds 1-6) in ISeta et al] (Il998h . The mass is 
assumed to be uniformly distributed within the radius of the 
molecular cloud complex. Rue = 50 pc. Given the spheri- 
cal symmetry of the model, we combine the 7-ray spectra of 
SRC-landSRC-2. 

Effectively, our model has three adjustable parameters, Z)28, 
^esc^MC, and s. As shown in Fig.|4l we determine /^esc^MC 
(i.e., the normalization) and s to reproduce the 7-ray spectrum 
for three values of the diffusion coefficient: D2S = 0.1, 1,3. 
For slow diffusion of D2S = 0.1, a steep index of ^ = 2.6 is re- 
quired since CR protons with p < 400 GeVc~^ are still within 
the GMC complex. CRs at higher energies are leaking from 
the complex, and correspondingly the 7-ray spectrum exhibits 
a break at ~ 50 GeV. Though the steep injection index re- 
quired for the slow di ffusion case is not theoretically expected 
(ICaprioli et al.l |2010|) . this problem may be alleviated by 
adopting a larger value of 6. The total amount of CRs released 
into interstellar space amounts to Wesc ^ 0.3 x 10^^ (Mmc/5 x 
10^Mq)~^ erg. For a nominal value of Z)28 = 1, we obtain 
s = 2.0 and Wesc ^ 1.1 x 10^^ (Mmc/5 x lO^Mo)-^ erg; the 
7-ray spectrum in the LAT range is steepened due to energy- 
dependent escape from the GMC complex. Finally, we find 
s = 2.0 and Wesc ^ 2.7 x 10^^ (Mmc/5 x 10^ Mo)-^ erg to re- 
produce the spectrum in the case of Z)28 = 3. 

We estimate that the energy channeled into the escaping 
CRs amounts to Wesc ~ (0.3-3) x 10^^ erg. On the other 
hand, the CR content trapped in W44 is estimated roughly 
as Wp '^ 0.4 X 10^^ erg assuming a uniform CR d istribution 
between the SNR shell and radio filaments ( ^4.11) . A com- 
bination of the Fermi-LAI data to be obtained with longer 
exposures, and future operation of the Cherenkov Telescope 
Array (CTA) will put more stringent constraints on the mod- 
els through determination of 7-ray spatial distributions as a 
function of energy. 
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